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Abstract 

We explicitly compute the semi-global symplectic invariants near the focus- 
focus point of the spherical pendulum. A modified Birkhoff normal form 
procedure is presented to compute the expansion of the Hamiltonian near the 
unstable equilibrium point in Eliasson- variables. Combining this with explicit 
formulas for the action we find the semi-global symplectic invariants near the 
focus- focus point introduced by Vu Ngoc [33j . We also show that the Birkhoff 
normal form is the inverse of a complete elliptic integral over a vanishing cycle. 
To our knowledge this is the first time that semi-global symplectic invariants 
near a focus-focus point have been computed explicitly. We close with some 
remarks about the pendulum, for which the invariants can be related to theta 
functions in a beautiful way. 

1 Introduction 

The spherical pendulum is a paradigm for an integrable system with a focus-focus 
point. The classical treatment of the spherical pendulum can be traced to the be- 
ginnings of mechanics in the work of e.g. Lagrange and WeierstraB, in particular 
because it appears as a Umiting case of the Lagrange top. Whittaker [Ml paragraph 
55(ii)] gives the classical treatment, including the explicit solution in terms of el- 
liptic functions. Its reappearance in a modern context is related to the discovery 
of Hamiltonian Mondromy by Cushman and Duistermat |T2] . This non-trivial 
global property of the actions and the rotation number of the spherical pendulum 
has not been noted in the classical treatment. The non-trivial topology of cer- 
tain torus bundles in the spherical pendulum is best understood with the help of 
singular reduction, and the detailed treatment of the spherical pendulum in this 
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formalism is given in [HI chapter IV]. Within the context of KAM theory the fre- 
quency map of the system was shown to be non-degenerate for all regular values of 
the energy- momentum map by Horozov [23] (also see [2T]). Later in [24J Horozov 
showed that the frequency ratio map does have singularities, so that the isoenergetic 
non-degeneracy condition of the KAM theorem fails on codimension one curves of 
non-critical values in the image of the energy-momentum map F . The Hamiltonian 
monodromy has also been studied from a complex analytic point of view in [5] and 

i- 

In view of these works the spherical pendulum is one of the best understood 
Liouville integrable systems. In this article we want to add to this knowledge the 
semi-global symplectic invariants of the focus-focus point of the spherical pendulum, 
and the expansion of its rotation number and period near this point in terms of the 
invariant and the Birkhoff normal form. 

Well known smooth invariants of a dynamical system are the eigenvalues at 
an equilibrium point, which are invariant under smooth coordinate transforma- 
tions. For Liouville integrable Hamiltonians the local symplectic classification of 
non-degenerate singularities was given by Eliasson in [20]. More global invariants 
have also been considered. There are Bolsinov's orbital invariants [HI E] which are 
invariant under smooth conjugacy. Using these it was shown that the flow of the 
Euler top is C° orbitally equivalent to the geodesic flow on some ellipsoid |8]. More- 
over, this equivalence cannot be made [7]. Related results on C" equivalence of 
divergence free vector fields are given in p6] . Restricting to smooth symplectic coor- 
dinate transformations the corresponding invariants (in one degree of freedom) were 
introduced by Dufour, Molino and Toulet [TB]. The semi-global symplectic invari- 
ants for the Liouville foliation of a Liouville integrable Hamiltonian system with two 
degrees of freedom in the neighbourhood of the separatrix of a simple focus-focus 
point were introduced by Vu Ngoc [33], and for hyperbolic-hyperbolic equilibria by 
Dullin and Vu Ngoc [19]. The latter paper also contains and explicit calculations of 
semi-global symplectic invariants near a hyperbolic-hyperbolic equilibrium. 

Consider two Liouville integrable systems on the manifold with symplectic 
structure Vl and with energy-momentum maps -Fi, -F2 : M'^ — )■ M^, each with a simple 
focus-focus point. If these integrable systems have the same semi-global symplectic 
invariants at the focus-focus point then there is a symplectic map ^ : — )■ in a 
full neighbourhood of the separatrix of the focus-focus point and a map $ : — t- 
reparameterizing the integrals such that 

Fi o ^ = $ o F2 . (1) 

The invariants are given as the coefficients of the power series of the analytic part of 
the singular action expressed as a function of certain regular actions. Coincidence 
of the invariants is a necessary condition that the systems can be globally equivalent 
in the sense of Q. The coincidence of these invariants is also a necessary condi- 
tion for the corresponding quantum systems to have the same eigenfunctions. The 
eigenvalues do not enter this description. The classical analogue of this statement 
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is that the Hamiltonian or the second integral can be changed with the function $ 
without changing the invariants. When considering a weaker notion of equivalence 
which does not allow the transformation $ then the additional invariants are the 
coefficients of the (simultaneous) Birkhoff normal form of F . 

The main objects in the classification are two different kinds of actions. First 
of all the standard canonical action with its singular behaviour near a separatrix 
is essential. But in addition also "imaginary actions" which are sometimes associ- 
ated with the quantum mechanical tunnelling play a crucial role. In the spherical 
pendulum they are given as complete elliptic integrals over cycles (3 and a that 
form a basis of cycles on the corresponding complex torus. It turns out that the 
imaginary actions are intimately related to the Birkhoff normal form near an un- 
stable equilibrium point, more precisely they are the Eliasson coordinates defined 
in a neighbourhood of the equilibrium. The key idea for the construction of the 
semi-global symplectic invariants is to express the regularised canonical action as 
a function of the imaginary action. In this way naturally a symplectic invariant is 
obtained that is independent of the Hamiltonian. 

In [18] the information contained in the transformation of the Hamiltonian into 
normal form [HT| [201 EHl [2Sl ES] combined with the semi-global symplectic invariants 
[53] where used to find the form of the frequency ratio map near a simple focus- 
focus point. Here we apply the general theory to the spherical pendulum, and we 
show how the two essential power series, the Hamiltonian in terms of Eliasson's 
coordinates and the semi-global symplectic invariant contained in the action, can 
actually be computed using a combination of adapted Lie-series normal form and 
classical analysis of elliptic integrals of the third kind. As a special case the semi- 
global symplectic invariants of the ordinary pendulum are found. We show that 
these invariants are given by the power series of Jacobi's nome in terms of the 
imaginary action. 

To our knowledge the present work gives the first example where the semi-global 
symplectic invariants near a focus-focus point are calculated explicitly for a partic- 
ular system. We hope that other systems will be added to the list, so that this will 
lead the way to a classification of Liouville integrable systems. 



2 Spherical Pendulum 

The spherical pendulum is a point of mass m constrained to move on a sphere under 
the infiuence of gravity. Let r G S*^ C be a point on the sphere ||r|| = I of 
radius I where the norm is induced by the standard Euclidean scalar product (■,■). 
Instead of the usual (singular) spherical coordinates on the sphere here we start with 
a global description that is free of coordinate singularities. With —rngSz being the 
force of gravity, Newton's equations of motion are 

mr = —mgez — Amr (2) 
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and the Lagrange multiplier A is 
r,r) -g{r,e^) 



A 



r, r 



determined such that d^(r,r)/dt^ = (r, r) + (r, r) = 0. The initial condition must 
satisfy d(r, r)/dt = 2(r, r) = in order to be contained in the tangent space of the 
sphere. The system has global Hamiltonian (shifted by a constant so that the upper 
equilibrium has energy zero) 

with momentum p = mr conjugate to r, and angular momentum L = r x p. This 
Hamiltonian with standard symplectic form dr Adp has (r, r) and (r, p) as constants 
of motion, such that the constraints (r, r) = l"^ and (r, p) = are automatically 
preserved. Using these constraints Hamilton's equations of motion become 

1 1 mgl mgl 

-L X r, p = — -L X p - — — + T— — (e^,r)r 



mP mP 

and are equivalent to The angular momentum = (LjC^) = xpy — yp^ is 
another constant of motion. Fixing the constrains defines a system on T*S'^ and 
this system is Liouville integrable because the Poisson bracket of the remaining 
integrals vanishes, {H,Lz} = 0, and H and Lz are independent almost everywhere. 
Formally the Hamiltonian ^ also describes the Lagrange top written in the space 
fixed frame and with the body symmetry reduced, see, e.g. [15]. The point r is a 
point on the symmetry axis of the body, and m/^ is the value of the equal moments 
of inertia. However, for the Lagrange top L is introduced as an independent variable 
(with the induced Poisson structure), and for the general motion (r,L) 7^ 0, which 
is impossible for the spherical pendulum. 

The frequency of small oscillations around the stable lower equilibrium ai z = —I 
is v = \fgjl- Measuring length in units of /, mass in units of m, and time in 
units of g all constants can be removed from the Hamiltonian ([s]). Since the 
transformation to non-dimensionalised variables is not symplectic we prefer to keep 
the constants for now. In the following we are interested in the dynamics near the 
unstable upper equilibrium at z = I. 



3 Linear normal form 

For the normal form calculation near the unstable equilibrium a local symplectic 
coordinate system on T*S'^ is needed. In the global Hamiltonian ([s]) we introduce 
local coordinates {x,y,l) in M'^ by {x,y,z) = {x,y, ^JP — x'^ — y"^) near z = I. The 
conjugate momenta {Px,Py,Pi) are found from the inverse transpose of the Jacobian 
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of the transformation of the coordinates and with pi = the obvious result that px 
and Py are unchanged and = is found. In these symplectic coordinates with the 
Hamiltonian reads 

H=^2 {pK^' - + pU^' - - ^^yP=^Py) + mg{^P-x^- y^ - I) . 

This Hamiltonian could have been derived directly from a Lagrangian using x and 
y as local coordinates. 

Before transforming the quadratic part of H into Williamson normal form [3] 
the symplectic scaling x = ^/y/mu, px = p^y/mu, similarly with r],Prj, is done. Thus 
we have proved 

Lemma 1. There exist local symplectic coordinates near the unstable equilibrium of 
the spherical pendulum with Q = dC, A dp^ + drj A dp^^ such that 

H = y(^^{pl+ pD - ^{^P^ + VPv? + ^( V^^^ - 1)) (4) 

where = + rj"^, = g/l and = mt^v = mgl/u . 

Taylor expansion of the potential (the last term in (|4])) gives 

_1.2 _ 1 4 _ 1 2 6 _ « f - (2n-3)!! . 

n=l 

SO that K, keeps track of the order of the non-linear terms. By scaling all variables 
^,r),p^,Prf with l/-\/ft (which is a symplectic transformation with multiplier 
all constants could be removed from the Hamiltonian up to the overall factor z///t, 
which could then be removed by a linear scaling of time. The two fundamental 
constants left in the Hamiltonian are l/u with unit of time and 1/k, with unit 
of angular momentum or action (i.e. kgm^/s). Without changing the symplectic 
structure and/or the time it is not possible to remove them. The usual approach to 
non-dimensionalise all variables (including time) is not used here until later because 
the semi-global symplectic invariants should allow comparison of different systems, 
and only in one of them these units could be chosen freely. 

So far point transformations were used to simplify H. To achieve the Williamson 
normal form [35] of the quadratic part of H it turns out that a symplectic rotation in 
phase space is necessary. The Hessian of the Hamiltonian has eigenvalues ±z/ with 
algebraic and geometric multiplicity two. Hence the equilibrium is a degenerate 
saddle-saddle point. Since there are no coupling terms each saddle can be treated 
separately in its canonical plane. A rotation by 7r/4 brings the hyperbola (p| — ^^)/2 
into the normal form piQi, similarly for {ri,pj^). Thus we have 

Lemma 2. The Williamson normal form at the unstable equilibrium of the spherical 
pendulum is achieved by the linear symplectic transformation 

= Qi- Pi, v^Ps = gi + Pi, V2r] = q2- P2, V2pr, = q2 + P2 ■ 



5 



The Hamiltonian in the new coordinates reads 



H = u (piqi +P2q2 - i^^{q^ -P^f + - a/1 - Kp2 + 




) 



K 



where = p\+ pi, q"^ = qf + ql and = p'^/2 + — {piqi + P2<i2)- 

The quadratic part of the potential has been absorbed in the quadratic normal 
form terms H2 = i^{piq2 +P2Q2)- The remaining terms of the potential are of order 
4 and higher. Together with the quartic term from the kinetic energy the quartic 
terms are 



The angular momentum in the new variables reads qiP2 — q2Pi, but it does not 
enter the quadratic Williamson normal form. This is why the quadratic normal 
form gives a degenerate saddle, while from the point of view of the foliation of the 
integrable system it is a focus-focus point. The reason for this is that any function 
of the integrals H and is again an integral, and thus has the same foliation. In 
particular the linear combination H+uLz for any u has an equilibrium of focus-focus 
type, with eigenvalues ±z/ ± iu. Incidentally this also shows that this equilibrium is 
a non-degenerate singularity in the sense of Eliasson. In fact we have achieved more 
than the Williamson normal form, because the second integral is also simultaneously 
normalised. See Appendix |A] for a few more remarks about this. 

4 Non-linear normal form 

From the work of [211 [201 EOl ES] we know that there is a convergent Birkhoff normal 
form near a non-degenerate equilibrium of a Liouville integrable system. We are 
now going to compute this normal form for the spherical pendulum. 

The equilibrium is of focus-focus type and thus we know that the Hamiltonian 
can be expressed as a function of two commuting integrals Ji and J2, where 

Ji = qm + q2P2, J2 = qm - q2Pi ■ (5) 

At quadratic order the Williamson normal form has already achieved this, H2 = 
vJi- But the higher order terms still contain other combinations of the variables 
qiiq2iPiiP2- Using non-linear near-identity transformations the Hamiltonian can be 
normalized to arbitrary order. For this we need 

Lemma 3 ( "Integral- Angle" coordinates). The integrals Ji and J2 as given by (|5| 
are momenta of a symplectic coordinate system on with the plane qi = q2 = 
removed. The canonically conjugate variables are 
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hq^ - p'f -^{p' + q^- 2(pigi + ^2^2))' • 



Hi 
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where only 62 is an angle. 



Proof. The construction mimics the usual action-angle variables near an elliptic- 
elliptic equilibrium point, however, here we are near a focus-focus point. Formally 
going into the complex plane we define C = qi + iq2, vr = pi + vp2 so that Ji — i J2 = 
and 9 1 + 162 = Then 

d(^i + 162) A d( Ji - i J2) = d In C A dC7f = dC A dvr 

and taking the real part of the previous formula and using 3?(d^ A d7r) = dqi A 
dpi + dq2 A dp2 gives the symplecticity of the transformation. The transformation 
is defined outside the set qi = q2 = 0. The inverse of the transformation is given by 

qi = cos 62 e^^ , Pi = (Ji cos 62 — J2 sin 62) e~^^ , 
q2 = sin 62 e"^ , p2 = {Ji sin 62 + J2 cos 62) e""^ , 

so that = e^^^ and p'^ = e'^^^{Jf + J|). □ 

Inserting this change of coordinates in to the Hamiltonian of lemma [2] a new 
Hamiltonian is obtained which is a function of Ji,6i, and J| only. 

It is important to keep in mind that at this stage Ji is only equal to the true 
Eliasson coordinates up to quadratic order. Normalising the Hamiltonian order 
by order each step of the normal form procedure introduces new variables Ji,6i, 
and then the tildes are dropped. In the limit of infinite order the Hamiltonian 
thus becomes a function of Ji and J2 only. In our particular case J2 is unchanged 
in the process, since H is independent of 62, in other words J2 already is a global 
integral. Note that the conjugate angle 62 is in general changed by the normalisation 
process. We are now going to describe the details of a practical implementation of 
this iterative procedure. 

The Hamiltonian is independent of 62, because this is the angle conjugate to the 
global integral = J2. The identity p^q'^ = Jf + J2 (using the notation for p and q 
of lemma [2]) allows to eliminate occurrences of p^ in H using p = J exp{—6i) where 

= -|- J|. Thus we find for the lowest order terms that are non-linear in Jj 

—H, = --{J^e-''^ - e^'^f - -{.Pe-^'^ + e^'^ - 2J^f 
vn 8 32 

= -1 J4 -401 lj^j2 -201 1^2 1^2 1^ 201 _ 1 40i 

32 8 16 ^ 16 ^ 8 32 

Now we seek a perturbation theory that removes the 9i dependence in this and all 
higher order terms. The classical approach with mixed variables generating functions 
is quite cumbersome at higher order; it is briefly described in Appendix [Bj The 
better alternative is the Lie series approach, in which the (near identity) canonical 
transformations are generated by the time e map of some flow. The Hamiltonian 
of this flow is the generating function of the transformation which only depends on 
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the old variables. The transformation is obtained by integrating this flow, which 
is conveniently done by computing iterated Poisson brackets in a kind of Baker- 
Campbell-Hausdorff formula. 

The abstract formalism is well understood, see e.g. [27]. A particular version of 
this type of perturbation theory has to specify a solution to the so called Lie-equation 
(or homological equation) 

{H2,Wi] = H.,-K, 

where Hi are the terms of order i of the given Hamiltonian, and Ki are the terms 
remaining in the transformed Hamiltonian, and Wi is the generator of the trans- 
formation, all at stage i of the iterated transformation process. The range TZi of 
the homological operator {-^2; "} is the set of terms that can be removed within the 
chosen class of transformations generated hy W^. A particular normal form consists 
of choosing a complement to the range Si, which is the set of unremovable terms at 
stage i. Hence Ki G Si and Hi — Ki & TZi. If the homological operator {H2, ■} is 
semisimple Si is the kernel of the homological operator. In the classical example of 
Birkhoff normal form near an elliptic equilibrium point Ki is the average of Hi over 
the flow of H2, while Ri is the oscillating part. The solvability condition for Wi is 
that Ki — Hi has zero average. 

We assume that the Hamiltonian is in Williamson normal form of lemma [2] and 
written in terms of the "integral-angle" coordinates of lemma [3} hence H2 = vJi 
for the spherical pendulum. A higher order term Hi is in normal form when it is a 
function of (Ji, J2) alone. This is similar to the usual Birkhoff normal form near an 
elliptic equilibrium. The essential difference is that here Ji is not a canonical action 
(it does not have a periodic flow, in other words, di is not an angle). Define a family 
of classes of functions Vi as functions of the form 

k=0 

where Qk is a. degree k homogeneous polynomial. The operator {H2,-} : Vi — )■ Vi 
simply gives Wi 1— )■ udWi/ddi. This operator is simple and splits into kernel and 
range. The kernel consists of functions Si that are independent of 9i. Because of 
the particular form of the class of functions Vi only the identity is in the kernel for 
odd i, while for even i the kernel consists of polynomials Qi/2- The range consists of 
functions TZi where each term does depend on 61, in other words Vi with Si removed. 
Using this in conjunction with theorem VII. B. 4 in [27] we can now prove 

Lemma 4. The Lie series normal form with respect to {H2, ■} for the spherical 
pendulum is unique and the normalised Hamiltonian is a function of Ji and J| 
only. 

Proof. We have already observed that this operator acting on Vi is simple. For 
uniqueness we need to check {Vi,Vj} C Vi+j and {Si,Sj} = 0. The latter is 
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true since Ji and J2 commute. The former holds because the commutator yields 
a function in Pj+j-i C Vi^j. The transformation given in lemma |3] combined with 
lemma [2] shows that all higher order terms Hi are in Vi. Finally we need to solve 
the Lie equation, and show that it has a solution in TZi. In our case the Lie equation 
reads 

Hj — Kj = V— — . 

Since the left hand side is in TZi (indefinite) integration with respect to 9i gives 
Wi G TZi. The normalised Hamiltonian is a function of J| only because the orig- 
inal Hamiltonian has the discrete symmetry J2 — > —J2 as well. Thus in fact the 
polynomials Qk are also even functions of J2, and this property is preserved in the 
normalisation process. □ 

Note that there are no small denominators in this process, since there is just one 
angle involved. Moser showed in [28j that even in the general non-integrable case 
a loxodromic equilibrium point has a convergent normal form. The above theorem 
shows how the Eliasson normal form of the spherical pendulum can be computed 
using Lie series. In the resulting normalised Hamiltonian all the higher order terms 
are in iS,, i.e. they depend on (Ji, J2) only, as desired. Thus we obtain 

Theorem 5. The normal form at the focus-focus point of the spherical pendulum is 
^HiJ,, J2) = Ji + ^{Jl + ?,Jl) - ^Wl + 9^2')+ (7) 

+ ^(5^1 + + 334) - + 410JiV2^ + 2714) + 0(J6) 

Proof. The proof is obtained by straightforward but lengthy calculation along the 
lines given in lemma |4| and the details are omitted. Here we just give the first steps 
of this computation to illustrate the process. There are no third order terms, so 
= = W3 = 0. Next the 9i independent terms in H4 can be read off from (j6]): 

— ^4= -4+-4 

Ku 16 ^ 16 ^ 
and integration of the remaining terms gives the generating function whose flow 
removes these terms: 



K UK J 128 16 16 128 

The terms of third order in J are the unremovable terms in Hq + ^{H4, □ 

Remarks: 1) Setting J2 = the normal form of the ordinary pendulum at its 
unstable equilibrium point is obtained. 2) It appears that by the non-symplectic 
scaling Jj — )■ 32Jj all coefficients of H{Ji, J2) become integers. 3) The normal form 
procedure works for any integrable system with a focus-focus point that has a global 
integral J2, so that the higher order terms are independent of ^2- 
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5 Action Integrals 



The normal form calculation of the previous section is local, and hence only valid 
in a neighbourhood of the equilibrium point. The local normal form actions are 
denoted by J,, while the global canonical actions are denoted by /j. In the present 
case the integral of angular momentum is J2 = -^2- The semi-global theory is valid 
in a neighbourhood of the pinched torus which is the separatrix of the equilibrium 
point. In order to find the semi-global canonical action Ii (which in this case is 
actually global up to monodromy) we introduce spherical coordinates (0, 9, 1) by the 
point transformation {x,y,z) = {I cos (p sin 6, 1 sin (p sin 6, 1 cos 6). Here the inclina- 
tion angle 6 is measured down from the north pole. The resulting Hamiltonian is 



H=lf^iy (pI+pI^^) +-(cos^-1) 
2 \ ^sm.B'^J K, 



where ^ = corresponds to the maximum of the potential which has energy zero. 
Here = = J2- This is the most compact form of the Hamiltonian, which 
unfortunately has a singularity at 6' = 0. The singularity can be avoided using 
singular reduction, see Cushman and Bates [111 chapter IV] . The final action integral 
is identical in both approaches. The action integral is obtained by integrating pedO 
over the real /3 cycle, hence 



1 f 2Kh , „ , K'^io 



nhih,j2) = — f W— -2(cos0-l)-^d^ 
Ztt If. V ly sm 



^/V2(/,-c + i)(i-C)-5l^ 



where the dimensionless length ( = cos 6 and the abbreviations Ji = kIi, J2 = 'tj2 
and h = Kh/u for the dimensionless actions and energy were introduced, n The 
non-trivial action Ji of the spherical pendulum is thus defined on the elliptic curve 

r,,, = {(C, w):w' = P(C)}, P(C) = 2(1 -e){h + l-0- Ji . 

Denote the roots of P by Q where — 1 < Co < Ci ^ 1 ^ C2- The elliptic integrals 
needed in the following have been given in terms of Legendre's standard integrals in 
e.g. [29]. The modulus and parameters for all of them are 

,2 Ci ~ Co Ci ~ Co 

K = — , n± 



C2-C0 ±l-Co 
With these definitions the action Ji can be written as 

= co{ciK{k) + C2E{k) + C3+H(n+, k) + C3_H(ra_, fc)) 



^ These tildes are dropped in the following! 
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■2 -2 

= V^^-Co) ' = 1 + ^ - C2 = C2 - Co, C3+ = C3- = The 

/3-cycle encircles the the interval [Co, Ci] along which nP' > 0. 

Lemma 6. The action of the spherical pendulum near the unstable equilibrium as a 
function of h, the energy difference to the critical value, and the angular momentum 
j2 has the expansion 

32 

2nh{h,j2) = 8- 2n\j2\ + J2 Arg(/i + ijs) + J2) In + 

+ h + + 3,1) - + 43,|ft= + 39,-) + , . . 

The coefficient in front of the logarithm ( the "imaginary action") is 

MK J2) = h-^{h' + 3jI) + :^h{h' + hjl) - ^{hh' + A2jlh' + 2ljt) + • • • • 

Proof. For the computation of the expansion of Ji the complete elliptic integrals 
need to be expanded in the singular limit that the modulus — )■ 1. We set h ^ he 
and j2 — ^ 32^ and expand in e. The roots of -P(C) in this limit are 

Co = -l + \3l + 0{e') 



so that the modulus becomes 
1 



and the parameters are 



n+ = l + ^(/i- V^' + il) +0(5^), n_ = -^ + 0(1/5). 



The fact that n_ behaves like does not cause problems because C3_ has a factor 
j|. For the expansion in the present so called 'circular' case the integrals of third 
kind need to be rewritten in terms of Heuman's Aq function, see e.g. [9]. Since 
there are two integrals of third kind some simplification is achieved when they are 
combined using the addition law for Aq, again see [9]. The result is 

h = CoCiK{k) + CoC2E{k) - -^Ao{<^, k), 
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where 



. / j1/2-Ci-Co ~ , . j1 J2 / (C2-i)(C2-Ci)" . 

09 = TT— arcsm ^ / , Ci = /i+l— G ; r \ ; ^ sm 09 . 

V 4(1 + C2) 2y 2(C2 + 1) ^ 

From this point onwards the expansion of Ji to high order is rather lengthy but 
straightforward using the formulae given in [H |9]. For a few more details of a 
similar calculation in a different system with a focus-focus point see |16j . 

□ 

The non-smooth term — 27r|j2| in ^ is present in order to make the numerical 
values coincide with the definition as a real integral, along with the choice of the 
principal argument Arg. In particular ^2) is an even function of j2, whose 
j2-derivative limits to 1 for h > and to 1/2 for h < 0. Since the spherical 
pendulum has monodromy, integer multiples of j2 can be added to Ji in order to 
make it locally smooth. When Ii is made smooth in this way it becomes globally 
multivalued. Making use of the freedom to add integer multiples of 72 the term 
could also be omitted. 

The function Ji (which multiplies the singular In-term in Ji) is also a complete 
elliptic integral over the same differential pgdd, but along the a-cycle that is van- 
ishing in the limit approaching the equilibrium. The a-cycle is a closed cycle in 
the original 6 variables. When the integral is rationalised hj z = cos 6 this cycle is 
folded once over itself, so that in the (^-variables it must be multiplied by 2. This 
factor of two is absorbed in our definition of a, which traverses the cycle twice in 
the ( variables. With this convention we have 

where a goes twice around the interval [Ci;C2] along which w'^ < 0. In the limit 
e — )■ these roots coalesce at 1. Since w is pure imaginary along this path we 
multiply by — i to make Ji real (by abuse of notation nevertheless we call Ji the 
imaginary action). 

As a real integral Ji is divergent, but as a complex integral it is well defined, 
and it can be transformed to Legendre normal form. The real expression in this so 
called 'hyperbolic' case can best be expressed in terms of Jacobi's Zeta function [Ij. 

Instead of using the complicated Legendre normal form, for Ji the series expan- 
sion can be found more easily from an expansion of the integrand. In contrast to the 
integral Ii an expansion works for Ji because the integration cycle is vanishing in 
the expansion limit. The residue of the differential form u = w / {1 — (^'^)d( at ( = ±1 
is =F|j2|/2. When h = j2 = the polynomial factors as P{() = 2{( + — 1)^ so 
that in this limit the pole of a; at C = 1 cancels. When h and j2 are both small the 
differential u can be expanded first into a Taylor series in e after setting h ^ he and 
/ — )■ le, and second in a Laurent series at ( = 1. Each term is polynomial in h,j2 
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and rational in (, y/( + 1 so that the residue at C = 1 of each term can be easily 
computed. Combining the two expansions an explicit formula for the coefficient of 
terms of degree n in h,j2 is found to be 

1 / (C - 1)" 9" w 

^(n - l)!9C"-i V d^l-C 

This gives the degree n term in the series of Ji(/i,j2) as stated in lemma [6] The 
factor of two appears because a is traversed twice in the C variables. 

In a semiclassical context Ji/2 may have the interpretation of a tunnelling in- 
tegral, which is used in uniform WKB quantisation across a separatrix. For a tun- 
nelling integral the integration path in the original variable 6 is not closed, but goes 
from one side of the energetically allowed region to the other. For our imaginary 
action Ji the cycle is closed in the original variable 6, thus the factor of 1/2 appears 
when considering the tunnelling integral. The interpretation of the imaginary action 
as a tunnelling integrals does hold for the ordinary pendulum, which essentially is 
obtained from the spherical pendulum by setting j2 = (see section [s]). For the 
spherical pendulum, or in general for a system with a focus-focus point, however, 
the uniform semiclassical quantisation is more complicated, see fIU{ |32] ■ In [2S] a 
semiclassical treatment of the spherical pendulum with an added cardanic frame was 
considered. In its uniform quantisation tunnelling integrals are used, however, they 
diverge in the limit of vanishing frame. Thus the analogy of the imaginary action 
(the integral over the a-cycle) with the tunnelling integral is rather metaphorical 
for the spherical pendulum. 

From the classical point of view the interesting observation is that Ji{h,j2) is 
the inverse of H{ji,j2) = h with respect to ji, so that H{Ji{h, j2), 32) = /i is an 
identity. Thus we have 

Theorem 7. The Birkhoff normal form near the focus-focus point of the spherical 
pendulum as given in Theorem is also given by the inverse of the elliptic integral 
Ji{h,j2) (over the a-cycle which is vanishing near the focus-focus equilibrium point) 
with respect to its first argument. 

Proof. After reduction, the spherical pendulum with one degree of freedom 9 is given 
by ([s]) . The dependence on the parameter j2 is ignored in the following. Recall that 
the transformation to (ordinary) action angle variables near an elliptic equilibrium 
point in one degree of freedom can be found from the mixed variables generating 
function S{I,q) = f^pdq, where p is obtained by solving H{p,q) = h{I) for q. 
This generating function is fixed by the requirement ^ d0 = 2tt, and this gives 
f pdq/{2iT) = I where the closed loop integrals are performed over the real solution 
of the equation H{q,p) = h (see, e.g. [3]). Thus the series expansion of the inverse 
of I{h) gives the Birkhoff normal form at the elliptic equilibrium point. 

The only modification we have to make to adapt this construction to our case of a 
hyperbolic equilibrium point is to allow a closed loop integral over a complex solution 
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of the equation H{q,p) = h. In particular consider a Hamiltonian of the form p^/2 + 
V{q) and consider a solution for which p is purely imaginary, i.e. the real integration 
range in q is an interval where no classical motion is allowed. The corresponding 
normalisation then is ^ d0 = 27ri, and accordingly ^pdg/(27ri) = J. The variable 
that is conjugate to J is defined by = dS/dJ as before. Performing this 
construction at the quadratic hyperbolic equilibrium H = \{a?'p'^ — h'^q^) reproduces 
the transformation to (Ji, ^i) as described in lemma |3} In the general case we have 
constructed Ji{h), whose inverse gives the Birkhoff normal form at the hyperbolic 
equilibrium point. 

Finally there is nothing special about a closed non-contractible loop that has real 
q and imaginary p; in the complex plane we can deform this path without changing 
J. In particular when the energy is above the energy of the saddle point of the 
Hamiltonian a path that encircles critical points with imaginary q and real p also 
works. □ 

The main advantage of using this theorem for the computation of the Birkhoff 
normal form instead of the Lie-series approach is that it is much simpler to first 
compute Ji(/i,j2) from the local expansion of the corresponding integral using the 



residue expansion (10), and then find the Birkhoff normal form by inversion of the 
series. 



6 Semi-global Symplectic Invariants 

From the general theory [33] it is known that the non-trivial action near a non- 
degenerate focus-focus point has the following form 

2vr/i(ji, j2) = 27r/io - In j - j) + S{juj2) 

where /lo is a constant and j = ji + ij2- Choosing a different branch of the complex 
logarithm changes Ji by an integer multiple of j2, and hence is a manifestation of 
the Hamiltonian monodromy of the focus- focus point. In the following theorem two 
different choices of the branch cut are made depending on the sign of j2, so that the 
result is consistent with 

Theorem 8. The non-trivial canonical action Ii of the sphercial pendulum near the 
unstable focus-focus point as a function of ji,j2 is given by 

2vr/i(ji,j2) = 8 - 27r|j2| + J2 Arg j - ji In |j| + ji + 5(ji,j2) (11) 
where the semi-global symplectic invariant S{ji,j2) is given by 

S{J,,J2)=J,\n32+^iJl+3J|)-^J^i5f,+5lJi) + -^ 
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Proof. We have already noticed in Theorem [7] that Ji{h,j2) is the inverse of the 
Birkhoff normal form, and thus inserting the normal form ^ into the series expan- 
sion of the action ^ produces the term —ji In \j \ , as desired. The miracle in this 
computation is that the coefficients of the regular part of j^) as obtained from 
the expansion of the elhptic integral are rational functions of h and j2, see (|9]), but 
upon substitution of the Birkhoff normal form ([T]) they turn into polynomials in ji 
and j2- The general theory [33] predicts this result, but it seems to be hard to see 
this otherwise. The resulting regular power-series part 5'(ji, 72) of the series /(ji, j^) 
is the semi-global symplectic invariant of the spherical pendulum at its focus-focus 
point. □ 

Remarks: 1) Setting j2 = the semi-global symplectic invariants of the ordinary 
pendulum at its unstable equilibrium point are obtained. 2) As noted in lemma |4| H 
is a function of j| only, and the same action Ji is also even in j2. Thus the discrete 
symmetry of the Hamiltonian becomes a discrete symmetry of 5*. 3) The coefficients 
of S (except for In 32) appear to become integers when each momentum is scaled 
by 32. The fact that there are rational coefficients appears to be related to the fact 
that the spherical pendulum is analytic. 4) The maximal absolute error between 
the action integral Ii{h,j2) and the approximation /i(ji,j2) given in Theorem [s] 
together with the formula for Ji(/i, 72) of lemma |6] using only the terms displayed in 
these formulas is less than 0.0001 inside a circle of radius 1/2 centered at the critical 
point {h,l) = (0,0), and less than 0.0032 inside a circle of radius 1. 5) In this and 
the previous section we are using scaled quantities which absorb the constants u and 
K. Unsealing shows that u does not appear in the symplectic invariant, while 1/k 
with units of an action naturally does appear in order to give 5* the correct units 
of an action. But the unsealing also gives the term jilnj>c in Ji, and one could 
argue that expanding the In gives another contribution to S as ji In k. However, the 
argument of In should be dimensionless, and thus one should leave the k inside, and 
do not incorporate it into the symplectic invariant. 

In [18] it was shown that the rotation number has the form 

27rW{ji,j2) = - arg j - A\n\j\ + ASi - S2 mod 27r , 

where Si = dS/dji. The function A can be computed from the Birkhoff normal form 
of H given in ([T]) as A = dj^H/dj^^H. Notice that for an elliptic equilibrium point 
where both ji are true actions the function A would be the rotation number itself. 
In the present case it is a major ingredient of the rotation number, but the rotation 
number also depends on the foliation via the semi-global symplectic invariant 5", and 
it has leading order singular terms which are universal for focus-focus points. After 
expansion of the quotient the result for A is 

3 15 15 45 

Aiji,j2) = -j2-Y^Jij2 + Y^j2(3j? + 2j2)-^^jij2(25jt + 43j2) + ... (12) 
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In order to obtain the rotation number given by the complete eUiptic integral ob- 
tained from —dli{h,j2)/dj2 the necessary choice of branch cut gives 

27rM/(ji, J2) = -^^li^ = 27rsgni2 - Arg j - A\n \j\ + AS, - . 
0J2 

The additional term preserves the property that W is odd in j2 and that the limiting 
values on the axis j2 = are ±1 for ji > and ±1/2 for ji < 0, since for = we 
have A = S2 = 0. This gives a very compact formula for the rotation number of the 
spherical pendulum near the focus-focus point. |^ 

The negative rotation number is obtained from the action by partial differen- 
tiation with respect to j2 keeping h = const. Differentiating ^ we see that the 
function A multiplying the In-term in the rotation number is dJi{h,j2)/dj2- This 
shows that A is the integral over the same differential as W, but instead over the 
vanishing cycle. Hence A is a complete elliptic integral as well. Moreover, this shows 
that the series of A{h,j2) can be found by differentiating the series of Ji(/i, 72) with 
respect to j2- Another even simpler pair of functions related to (Ji, Ji) are the 
periods of motion obtained by differentiating {Ii{h,j2), ^2)) with respect to h, 
instead of with respect to 72 • In all cases the series expansions of the integrals over 
the vanishing cycle can be easily found by expanding and computing residues as 
described for Ji. 



7 Dynamical Interpretation 

The physical interpretation of the rotation number is easily understood when con- 
sidering a rational rotation number W = p/q. This means that the orbit completes 
q periods in the inclination angle 6, while it completes p periods of the azimuthal 
angle (p around the pole. 

From a dynamical point of view the most interesting quantities that have simple 
interpretations are derivatives of the action Ii{h,j2) which give the period of the 
(symmetry reduced) motion and the rotation number. The period is found to be 

. _ dh/dn _ - In |j| + ^1 _ 1^ § + il^i + 

where Si = dS/dji. The period diverges logarithmically when approaching the 
unstable equilibrium point at j = 0. The ji derivative of the symplectic invariant S 
gives the smooth correction to this singular behaviour, when the factor l/{dH/dji) 
multiplying the singular term ln|j| is factored out. The original construction in 
[33] achieves the splitting into singular and smooth contribution using two Poincare 
sections transverse to the stable respectively unstable manifold of the equilibrium. 

^ In Cushman and Bates the rotation number W is an even function of j2 , which is obtained 
from our odd W by taking the absolute value. Thus W becomes continuous and odd. 
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Figure 1: Periodic orbits of the spherical pendulum in stereographic projec- 
tion of the sphere fixing the unstable equilibrium point with initial condi- 
tions along the circle in (ji,j2) space with radius 0.75 with rotation numbers 
4/7, 3/5, 5/8, 2/3, 5/7, 3/4, 4/5, 7/8. Individual pictures are shown at the same scale 
cropped to [—8, 8]^. 

This divides the separatrix into an inner and an outer part, where the inner part 
approaches the equilibrium point for \j\ — )■ 0. This gives an interpretation to the 
coefficient In 32 of ji in S, which is the leading order term of the time it takes to 
traverse the outer part of the separatrix. 

In the plane of (j2;Ji) "we now introduce polar coordinates by j2 = rcoss, 
ji = rsins. The roots become Co = ~"1 + |r^cos^s + O(r^) and Ci,2 = 1 + 
± sin s)(±r)(l ± ^r) + 0{r^). Recalling that ( 

= cos 6' shows that the min- 
imal distance of the orbit to the north pole measured in the inclination angle 6 is 
\/r{l — sin s) + 0(r^/^), while the minimal distance to the south pole is r| coss|/2 + 
O(r^). For s — > 7r/2 the angular momentum j2 vanishes and both distances vanish. 
For s —7- — 7r/2 the angular momentum also vanishes, but since ji ~ /i < the north 
pole is still inaccessible. In figure [T] the radius of the excluded region in the centre 
of each picture is \/r{l — sins)/2 to leading order (the factor 1/2 comes from the 
stereographic projection), while the overall size is 4/(r| coss|) to leading order. 

It is interesting to note that to leading order the modulus fc^ = 1 — r/2+T^(4 + 
3sins)r^ + O(r^) is independent of s. The leading order of the rotation number W 
for J2 > in polar coordinates is 

3 s 3rcoss , , 32 , 

This clearly shows that along the half-circle s G [— 7r/2,7r/2] in the right half-plane 
with j2 > the rotation number changes from W = 1 /2 to W = 1, so hj continuity 
every rotation number in between occurs at least once on each such half-circle. At 
least for small r every rotation number occurs exactly once along the half-circle with 
radius r. 

Another quantity that is also interesting is the twist T given by the derivative of 
the rotation number with respect to j2 for constant h. The zeroes of this function 
indicate vanishing twist, see [18]. Using the general formula derived there for the 
spherical pendulum we find 

ii Q 3 32 
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Figure 2: Two distinct periodic orbits of the spherical pendulum with same energy 
h = 0.05 and same rotation number W = 5/6 shown in stereographic projection of 
the sphere fixing the unstable equilibrium point in three different scales. Existence 
of two orbit with the same rotation number for fixed energy implies that the twist 
vanishes for some initial condition "in between" . 

where Wi = dW{ji,j2)/dji. In polar coordinates the twist becomes 

^ sin s 3 , 32 15 3 ^, , 

2nr = + -ln — --cos2s + Or . 

r 8 r 16 8 

This shows that for r — )■ to leading order the twist vanishes along the line s = 
0, along which the rotation number approaches 3/4. Considering the next order 
correction one can see that for ji > only rotation numbers in the interval (3/4, 1] 
occur. The equation T = can be approximately solved for s. This shows that the 
curve of twistless tori curves upward from the horizontal axis ji = 0. The parametric 
form in polar coordinates shows that at the origin the twistless curve has slope zero, 
and that higher order derivatives do not exist. 

Evaluating the rotation number on the curve of vanishing twist and expanding 
for small r gives 

,,,, 3 3r, 32 5, 

W* ^ - ^ (In ) . 

4 87r^ r 2^ 

8 The Pendulum 

The pendulum Hamiltonian is obtained from the spherical pendulum Hamiltonian 
(§ by setting j2 = 0:f\ 

1 2 ^ / n \1 2 ^ 2 ^ 

H = -Kupg + - (cost' - 1) = -Kupg - 2- sin - . 

•^For the spherical pendulum we have chosen h — at the unstable equilibrium point which 
corresponds to 9 — 0. For the pendulum other conventions may seem more natural, but we stick 
with this unusual choice for ease of comparison with the spherical pendulum. 
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Introducing the scaled energy and action as before the canonical action / = Ii{h,j2 = 
0), the imaginary action J = Ji{h,j2 = 0) (which does have the interpretation of 
twice the tunnelling integral in this case), and their /i-derivatives T and U, the real 
and imaginary period, respectively, for /i > are given by 

I^(h) = — (f \/2{h + 2sm^9/2)d9 = -^E{k), A;^ 



27r V ^ ' ' nk ' 2 + h 

J^{h) = — i \l2{h-2sin\i^e/2)de = \{K{k') - E{k')) 
271 J a V Ilk 

TJh) = 27r^^ = 2kKik) 
oh 

U^^h)=27i^^ = 4kK{k'), k'^ = l-e- ^ 



dh ' " 2 + h 

For h < the corresponding results may look unusual because they are not the true 
action or period of the pendulum, but the analogue of the action and period for 
the spherical pendulum, which differ by a factor 1/2 (and thus the action becomes 
continuous across h = 0). The results are 

i4h) = -{E{k)-{i-e)K{k)), e = ^^ 

J_(/i) = -2/_(-2 -h) = -{k^K{k') - E{k')) 

71 

T^{h) = 2K{k) 

U-{h) = 2T_(-2 - h)/27i = AK{k'), k'^ = ~>0 

The simple relation between 2T_ and t/_ given by swapping k for k' was first noted 
by Appell 0|. His interpretation of the imaginary period U is that it is obtained 
by replacing time t in the equations of motion by it. We see that a similar relation 
holds for the action 21 and the imaginary action J, except for another overall minus 
sign which is chosen such that j = h + 0{h^) for small /i, as required by the Birkhoff 
normal form near the hyperbolic equilibrium. It should be noted, however, that 
Appel's relations are only true for energies below the critical energy h = 0. The 
general question of the physical meaning of the action or period integrals over the 
a and (3 cycles in post quantum mechanics times is answered by uniform WKB 
quantisation, in which the action and the imaginary action both appear. In this 
paper we add another interpretation which is the relation between the Birkhoff 
normal form at an unstable equilibrium point and action integrals of the a-cycles 
that vanish at this equilibrium point. 



19 



The series expansions for both signs of h are 

2nl = 8 + h+ (^h - + . . . ^ ln(32/|/i|) + + ... 

J = h-—h^ + ... 
16 

T=i^l-h-^h' + ...yn{32/\h\) + \h + ... 

U/27r = l--h+ —h^ + ... 
' 8 256 

Remarks: 1) The four integrals satisfy the Legendre relation lU — JT = 8. 
2) Another way to interpret the actions given above for the case that < is to 
consider the pendulum reduced by the discrete symmetry {9 — vr, p) — )• {—9 + n, —p). 
Alternatively certain integer factors need to be introduced to recover the true semi- 
global symplectic invariants of the pendulum. See [19] for a discussion of these 
integers in the case of the C. Neumann system. 3) The two branches of J and U for 
h < and h > join analytically at h = 0. 

For the pendulum the calculations are simpler because they do not involve elliptic 
integrals of the third kind. In fact J{h) and I{h) are of the second kind and are 
simply the integrals over the a- and /3-cycle of the same elliptic curve. The semi- 
global symplectic invariants are obtained by computing Jo J~^, i.e. the integral over 
the /3-cycle as a function of the integral over the a-cycle. The result is 

27r/(j) = 8 + j(l + ln(32/|j|)) + + • • • 

This can be done efficiently to high order because the series expansions of the com- 
plete elliptic integrals K and E are well known. Hence we find the quadratic and 
higher order terms of the semi-global symplectic invariant of the pendulum are given 
by 

3f 5f , 55j^ I89f ^ 3689f _ 3129f ^ 1575405jS 



32 512 32768 524288 41943040 134217728 240518168576 
This can be obtained from the semi-global symplectic invariant of the spherical 
pendulum Theorem |8] by setting j2 = 0. 

We are now going to connect the semi-global symplectic invariant to the Riemann 
matrix and i} functions. The main observation is that if we differentiate I{J) with 
respect to J we obtain the ratio /'/ J' where the prime denotes differentiation with 
respect to the energy. Now /' and J' are the real and complex period, so that this 
gives the period ratio of the curve (up to a factor of i). |^ 

dl 31/ dh T K{k) ir' 1 32 ^ iJ + Oif\ , ^ ^ 



dJ dJ/dh U 2K{k') 2 2tx \ \j\ l-b' + 0(j2 



^When we speak of period ratio in this section we mean period ratio t of the ehiptic curve 
instead of the rotation number W . 
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where r' denotes the reciprocal of the usual period ratio r = iK'/K, since we are 
expanding in h, hence in k'"^ instead of in fc^. 

Notice that by passing from I{J) to its derivative we do not lose any information 
since the constant term is not part of the semi-global symplectic invariant. 

In more degrees of freedom and near completely hyperbolic points this may give 
a variant of the Riemann matrix, e.g. in the Neumann system [19]. The pre-potential 
computed in [19] is a function of imaginary actions whose Hessian is the Riemann 
matrix. In particular in this higher dimensional situation it may be advantageous 
to work with the Riemann matrix instead of with I (J) as invariants). 

Instead of the period ratio r (or the Riemann matrix) commonly the so called 
nome q = e^'^''^ is used. Passing to the nome converts the logarithmic singularity in 
the period T into a pole. Exponentiation thus leads to a power series in j 

q = exp — ZTT — 



dJj 

= l-6l^ + 48/2 - 436/^ + 4254/^ - 43452/^ + 458192/^ ± . . . (13) 

where I = j/32. Alternatively we may consider the reciprocal which still 'remembers' 
that there was a singularity which is now turned into a pole 

exp I 27r-^ ) = ^ + 6 - 12/ + 76/^ - 606/^ + 5412/^ ± . . . 
\ oJ J I 

The coefficients of either series are an equivalent set of "modified" semi-global sym- 
plectic invariants. Hence it appears that the nome itself (or the Riemann matrix in 
general) gives an equivalent set of modified semi-global symplectic invariants, when 
expressed as a function of the imaginary action J. This is a very natural object to 
consider from the complex analytic point of view. The expansion of the nome in 
terms of the modulus k"^ is well known (the inverse of which can be expressed in 
terms of theta functions). Now k'^ = 2/{h + 2) for the pendulum and the Birkhoff 
normal form allows us to replace h by j and thus gives the final modified semi- 
global symplectic invariants. Thus we can compute q{J) from known functions and 
the Birkhoff normal form (whose functional inverse is also expressed in terms of 
known functions). 

It turns out that even though we cannot find q{J) directly, its inverse func- 
tion J{q) can be expressed in terms of extremely fast converging theta series = 
^(— l)"g" . By standard identities involving thetanullwerte [Ij we find that 

J[q) = -16q ,3 = 8g— -2 = 8q 



^dq^l ^dq2K{q) 
q-Aq^ + 9g9 - 16gi6 + 0(^25) 
1 - 2g + 2g4 - 2g9 + 2q^^ + 0(g25))3 
32(g + + 24g^ + 76g^ + . . . ) 



32 



^The sign in the exponent is fixed by convention so that 3(t) > which imphes that the 
Riemann matrix has positive definite real part, so that q = = e"^^ — e^^^ 1^ < 1. 
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and the inverse of this series gives the modified semi-global symplectic invariants. 
The coefficients of the g-series of ti/{2K) and of J/32 have been studied by J.W.L. 
Glaisher [22] ■ He did, however, not consider the inverse of this series. Since the 
formal inversion of a power series is an invertible process, one may consider taking 
J{q) itself as the invariant instead of its inverse. However, this seems to far removed 
from the original construction, in which J provides a local coordinate system near 
the singularity. 

The above formula for J(g) incidentally also proves that the modified semi-global 
symplectic invariants are integers after scaling by 32, since division by a series with 
constant coefficient 1 and integer coefficients gives a series of the same form, and 
similarly for inversion of a series with linear coefficient 1. The fact that these 
coefficients are integers may be taken as another indicator why it may be a good 
idea to consider the series exp{—27idI/dJ) instead of the analytic part of /(J) as 
the semi-global symplectic invariants. 

Let us finally comment on how similar expressions can be obtained for the spher- 
ical pendulum. The key idea is to pass to the derivatives of the invariant, namely 
Ti = dli/dji and T2 = dli/dj2, which are entries of the period lattice, see [33|. 
Now observe that the singular terms are ti = — In |j| + . . . and T2 = arg j + . . . . 
Thus exp(— 27rri) = |j| -|- . . . and exp(27rir2) = j/\j\ + ■ ■ ■ so that the combination 
q = exp(— 27rri + 27rir2) is free of singularities and is a complex analogue of the nome 
q in the case of the pendulum. Repeating this construction with the full expression 
for Ji in ([TT]) including the symplectic invariant gives 



q = jexp{-Si + 15*2) 
= /(I + 6(/ - 21) - 3{17f + 211- 35F) + 2(37/^ + 666/^^- 633//"- 288/^) 



where / = j/32 = (ji + ij2)/32. For / = / this reduces back to (13) for the ordinary 
pendulum. The construction can be done with —Si— 182 as well. The main question 
that remains is whether the inverse of this series giving ji + ij2 in terms of q and 
its conjugate has a simple expression in terms of say '(?-f unctions. This is harder to 
do than for the pendulum because the integral Ji is of third kind. But expressions 
for complete elliptic integrals of the third kind in terms of the nome do not seem to 
be readily available, and would also need to depend on a second parameter, which 
would presumably be the imaginary part of q). 
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A Simultaneous linear normal form 



The symplectic linear normal forms near equilibria are due to Wiliamson [35J. Let 
the original variables be {^,p^,r],pr,) with symplectic form f2 = A dp^ + dr] A 
dpjj, and assume that the equilibrium point is at the origin. The Hessian of the 
Hamiltonian at the origin is denoted by D^H. Let J4 be the standard symplectic 
matrix corresponding to Q. For a linear focus-focus point the eigenvalues of J^^D^H 
are A = ±a ± iu, where we assume au 7^ 0,0; > 0. The stable eigenspace is 
spanned by the eigenvectors Vg^Vg of eigenvalues A^, A^ = « ± and the unstable 
eigenspace is spanned by Vu of eigenvalues A^, A„ = — a =F icu. In each eigenspace 
the dynamics is that of a (contracting or expanding) focus point, hence the name. 

According to the Williamson (i.e. linear) classification the spherical pendulum 
has a saddle-saddle equilibrium with degenerate real eigenvalue and trivial Jordan 
blocks. From the non-linear point of view of the foliation of the integrable system it 
is a focus-focus point, because we may form a linear combination of the constants of 
motion H and without changing the foliation. This is why we treat the saddle- 
saddle degenerate focus-focus case, as is appropriate from the non-linear 
point of view. We know that for the non-linear system the angular momentum 
is a constant of motion, even though it does not appear in the quadratic part. 

In fact we want a little bit more than the Williamson normal form, since we 
want to normalise the Hamiltonian and the second integral simultaneously. Since 
they commute, they can be simultaneously normalised. In the case of the spherical 
pendulum the angular momentum is already in normal form = C,Pr] — VPs.- Thus 
we seek a linear symplectic transformation M from (^,pg, r],prj) to {qi,Pi, ^'2,^2) that 
leaves J2 = C,Pr] — VPi = Q1P2 — (I2P1 invariant and transforms the Hamiltonian to 
aJi + UJJ2 where Ji = qipi + 52^2- The eigenvalues of the equilibrium at the origin 
are A = ±a; ± iu, which in our case simply gives A = ±z^, since w = 0. Hence M 
must satisfy 

M* J4M = J4, M^D^J2M = ^2 J2, M^D^HM = aD^Ji + ujD^J^ . 

Note that since {^1,^2} = also their Hessians D^Ji and D^J2 commute, and hence 
they diagonalise in the same basis. The transformation given in lemma |2] satisfies 
these conditions. 



B Canonical Perturbation Theory 

Instead of Lie-series perturbation theory the analogue of canonical perturbation 
theory may be used, which we briefly describe in this appendix. Considering as 
given in (|6]) as a series of exponential function we now introduce a formal averaging 
operator. Given a function 
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we define 



{F)^fo, {F} = F-{F). 

The "oscillating" part {F} can be integrated to give the generating function of 
a canonical transformation. The average (F) cannot be removed by a canonical 
transformation. But it is independent of 9i, and this is the desired outcome. The 
procedure is analogous to ordinary perturbation theory in which 9i actually is an 
angle, and the series is trigonometric instead. The formal analogy works on the level 
of the series expansion, while the usual interpretation as "averaging over the flow 
of the unperturbed system" does not work directly in our case because this flow is 
unbounded. This is so because Ji is not an action, i.e. it does not have a periodic 
flow. Nevertheless, Ji is a local integral near the equilibrium point, thus the name 
"integral-angle" variables. 

In the present case this procedure gives 

and by integration 

128 16 16 128 

The mixed variables generating function of the canonical transformation is 

S{Ji, J2, Oi, 62) — J\9i + J2^2 + £'S'i( Ji, J2, 9i, 62) 

where Ji are the new momenta. Inverting the transformations explicitly up to first 
order gives 

Ji^Ji + edSi{h,J2AM/dOi 

ei^9i-edSi{Jij2:0i:02)/dJi 

This leaves the angular momentum J2 unchanged, while Ji is modified by —e{H4}. 
Both angles are transformed in a non-trivial way. 

C Rotation Number 

The rotation number can also be computed directly from elliptic integrals, instead of 
from the normal form and the semi-global symplectic invariants. For completeness, 
and to show the consistency of the approaches here we give the expansion of the 
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rotation number in h and j^- We start with the expression in Legendre normal form, 
which is 



% 27,Jg(l-C-)w ^ ,rV2(C2 - Co) V 1 - Co 1 + Co 
Similarly, the period of the reduced motion is 

T{h,j,) = 271 = / - dC = -4^Kik) . 



dh Jp w VC2 - Co 
Expanding the elliptic integrals in the limit h, j2 gives 

Lemma 9. The rotation number of the spherical pendulum near the focus-focus 
point has the expansion 

3 / 5 35 \ 32 

27tW = 27T sgn j2 - Aig{h + ij^) + -J2 [l - —h + —{K' + jI) + . . . \ \n — 

J2 5h' + Qjl J2h 77h^ + 174jlh^ + 93jl 
8 ^2 256 

where h is the difference in energy to the critical value, j2 is the angular momentum, 
and (? = il + . 

The smooth function multiplying In is given by the residue series of the integrand 
of W, in the same way as Ji is the residue series of the integrand of Ji. Taking 
this function, substituting ([T]), and expanding in ji and j2 reproduces the function 
MjiJ2) given by Q above. 
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